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Summary 


The  concapt  of  coharanca  at  a  maaaun  of  tha 
linear  relationehip  between  two  or  mora  tima  tar  tat 
it  discussed.  Tha  definition  of  coharanca  in  larma  of 
tha  complex  croaa  power  apectral  density  matrix  it 
given  and  tha  relationship  between  pairwise  and  multi* 
pla  coharanca  and  channel  aignal  to  noita  ratios  it 
discussed.  A  sample  autistic  for  coherence  it  given 
while  the  details  describing  the  performance  which 
car  be  obtained  from  this  statistic  are  contained  in 
a  separate  report  by  these  authors1. 


1.  Introduction 

A  problem  of  interest  in  many  different  disc* 
iplines  is  that  of  determining  if  there  is  a  measurable 
relationship  (physical  causality)  between  two  or  more 
time  series.  In  addition,  one  would  often  like  to  ob¬ 
tain  a  quantitative  meaningful  measure  of  the  degree 
of  that  relationehip.  This  paper  describes  one  possi¬ 
ble  measure  of  such  a  relationehip,  the  coherence. 

The  most  common  measure  of  such  a  relation* 

•  hip  is  the  pairwise  or  multiple  correlation  coefficient. 
The  nature  of  the  correlation  coefficient  is  well  docu* 
mented  and  will  not  be  discussed  here  other  than  to 
note  that  it  is  not  a  function  of  frequency  and  may  be 
effected  by  linear  transformations  of  either  of  the 
me  senes. 

Tne  coherence  function  (magnitude- squared 
multiple  or  pairwise  coherence  function)  is  defined  as 
a  frequency-gcnendent  quantity  that  ranges  between 
aero  and  one4*  This  coherence  function  is  sero  if 
the  two  or  more  ergodic  time  eeriea  are  independent 
(uncorrelated,  if  Gaoaeian)  and  equal  to  one  at  any 
frequency  where  there  it  a  linear  transformation 
between  the  one  or  more  input  time  series  and  the 
output  or  refer snee  time  eenee. 


*<t)  »  I  -  r)df  ♦  v(t).  (n  x  1)  (1. 1) 

Note  that  becauee  of  physical  causality  require¬ 
ments  t(r)  is  *aro  for  all  y  lest  than  aero.  In  fact*  It 
will  be  sero  for  all  y  lose  than  some  positive  time 
which  is  the  time  It  takes  a  signal  to  travel  from  the 
•ource  to  a  sensor. 

The  true  value  of  the  coherence  between  time 
series  is  generally  an  unknown  quantity.  In  fact,  any 
maasure  of  the  relationship  between  two  or  more  time 
series  generally  must  be  based  on  time  traces  of  those 
eerie#.  The  functional  relationehip  between  the  time 
series  and  the  measure  or  estimate  of  coherence  is  a 
•ample  statistic  for  coherence.  The  assumption  that 
any  one  infinite  length  sample  of  each  eeriea  will  be 
enough  to  allow  us  to  estimate  the  coherence  exactly 
is  made  implicitly,  and  all  time  series  sre  assumed 
to  be  stationery  and  ergodic.  Unfortunately,  in 
practice,  one  is  given  only  a  finite  amount  of  data 
from  each  of  the  time  series.  In  this  case  the  ssmplc 
statistic  is  a  random  variable  distributed  about  the 
"true  magnitude -squared  coherence"  The  density 
function  for  the  sample  statistic  in  this  report  is  given 
for  a  number  of  specific  values  of  degrees  of  freedom 
(N)  and  number  of  time  series  (Ml  in  reference  1, 

Based  on  these  density  functions,  receiver 
operation  characteristic  (ROC>  curves  have  been  de¬ 
veloped.  These  curves  define  the  probability  of 
detection  versus  probability  of  false  alarm  for  a 
signal  of  a  given  true  coherence.  Curves  of  probabil¬ 
ity  of  detection  versus  true  coherence  for  fixed  levels 
of  probability  of  false  alarm  have  «l*e  oeen  developed 
for  many  degrees  of  frtedJm  and  number  of  sensors 
up  to  IQ.  These  are  all  reported  in  detail  by  these 
authors  ir.  reference  1.  The  fundamental  question  of 
the  relationship  between  input  signal  to  noise  levels 
and  true  coherence  for  the  two  channel  anc  multi¬ 
channel  cases  of  Figure  1  is  discussed  in  this  paper. 


The  situation  of  interest  is  shown  in  Figure  1 
where  P,  indicates  the  noise  contaminating  the  signal 
u(: '  in  the  Ith  channel.  In  general,  each  transmission 
channel  is  composed  of  linear  and  nonlinear  parts 
(Figure  2).  The  sum  of  the  output  of  the  nonlinear 
system  (usually  a  small  part  of  tne  total  transmission), 
tne  measurement  noise,  and  the  background  noise  is 
t  rouped  into  the  effective  noise  term  Vi(t)  (Figure  3). 

We  are  interested  in  detecting  the  presence  of 
a  common  signal  u(()  in  two  or  more  channels,  The 
nput- output  relationship  indicated  in  Figure  1.3  can 
oe  written  as 


For  s  general  discussion  of  the  concept  of 
coherence,  the  reader  is  directed  to  reference  5. 

For  a  detailed  derivation  of  the  distribution  of  tne 
pairwise  and  multiple  coherence  statistic,  the  reader 
is  referred  to  reference  7.  The  densities  and  derived 
performance  curves  ir.  reference  1  are  particularly 
difficult  to  obtain  for  low  coherence  and  high  values  of 
N  (number  of  samples  of  the  time  senes),  and  based 
on  the  authors  knowledge  are  not  available  elsewhere 
sn  the  literature. 
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I.  Croaa-Pow «r  Spectral  Danaity 
Mjlriij  miH  MwjrtgU  gjjwrwrg 

Multiple  coherence  can  b«  moat  eeeiiy  defined 
in  terma  of  the  croaa  apaetral  denaity  matrix 
S^lwl,  where  the  elemente  of  thia  matrix  are  defined 

V“,*svj,“,*y“>  (2-*' 

The  croaapower  apaetral  deneity  matrix  ia  of 
courae  equivalent  to  the  croaacorrelatioa  matrix. 

Cither  of  theee  together  with  the  meana  of  the  U 
jointly  gauaaiaa  etationary  proceaaea,  x.(t), x-(t), .., 

completely  apectfiee  the  joint  diafributton 
function  of  theee  proceaaea. 

Given  M  finite  length  time  tracaa,  there  are 
well  known  techniquee  far  obtaining  "temple  eatimatea" 
of  the  croaa-  and  autopower  apectral  elementa.  Theae 
aatimatea  are  uaed  to  obtain  aample  eatimatea  for  the 
multiple  coherence  between  the  varioua  time  aeriee. 

The  aample  eetimate  for  the  Croat -power 
apectral  denaity  matrix  ia  a  function  of  the  baaic  data 
x|  *2**''  ‘ ' '  ov*r  *em*  finite  time  record.  It 

doea  not  requireknowladf  «  of  any  of  the  character- 
iatica  of  the  tranamiaaion  channela  or  of  the  aignal-to- 
noiae  rattoa  of  the  received  elgnale.  The  meaning  of 
multiple  coherence  will  be  diacuaaed  in  terma  of  theae 
quantitiea  in  later  aactiona  in  an  attempt  to  Illuminate 
the  aubject.  Hare,  however,  we  will  define  multiple 
coherence  aimply  in  terma  of  the  croaapower  apectral 
cenaity  matrix  and  ita  elementa. 


Thta  it  the  mutual  or  pairwiac  coherence  between 
channote  one  and  two. 

In  the  three-channel  caee  it  ia  eeaily  aaan 
(writing  Sy(te)  ae  S(J)  that 

P'1,2.31  V  S  S  e  le  |Z  U“ 

SllS22S33'Sll|S2sI 


Uaing  the  fact  that  peirwlee  coherence  ie  defined 
by  equation  2.  $  we  can  write 


|n,2.jf  • 


h.2Nri./-2V^2523V/<SilS22S31> 

,  I  2  •' 


Ndte'that  if  the  croaapower  apectral  denaity  of 
channela  two  and  three  ie  aero  (S,3«0),  the  coherence 
between  theae  two  channela  ia  aero  and  the  multiple 
coherence  of  channel  one,  given  two  and  three  ia. 


I’i.of  •in.2|*+l>,i.jf 


The  multiple  coherence  between  x.(t)  and 
k|(t(,x-(t),...x.  ,(t),  x  ft). ... x  (t)  iaJ defined  by 
I  reference  S)  J*1  M 

|’j:1.2...j-l,J+l,...Hl2e  1-1/  [s..Ue)S-y(<e)l.  (2.2) 

unere  )  ia  the  jth  diagonal  element  of  the  inverae 

o:  ih«  Th«  multiple  coherence  of  the  jth 

•eneor  with  reepect  to  the  other  sensors  represents 
the  proportion  of  the  vartsnes  (power)  of  sensor  j  thst 
can  be  explained  by  a  linear  combination  of  tiut  remain* 
mi  sensors  in  a  minimum  mean  square  sense  . 

By  reducing  this  definition  to  the  simple  two- 
channel  case  we  can  write 


S22(W,'S12,W) 


S5*,“”*Su(u.)S2Z(w)-|Su(Uf 


-si2<«)Sn(«) 


1.  .  Pairwiaa  Coharanca  and 
Signal-to-Noite  Ratiot 

If  x]  (t)  and  x2(t)  ara  ganarattd  aa  indicated 
in  Figure  3.  equation  1 .  1  can  be  written  in  the 
frequency  domain  by  Fourier  treneform  ae 

XjlWlaCjlMtUflU.'aVjlW)  C 


XZ(W)=C2(U)U(W)+V2(U!).  (3.1 

Again  aa  turning  the  noiac^terma  Vj  and  V  are  inde¬ 
pendent,  the  croee-  a  lid  autopower1  apectral  dentitiea 
can  be  written  ea 

Su(ui).  Gj(iu)  2Su(io)+Svl(iu)  (3.  3 

S22(ut+ G2(u)  2Su(ei)+Sv2(ui)  (3.4 

|S12("l|i*  °1,W!2  <3.  S 


S1:(u,)*S22(o,).'[Sn(W)S22(«)  -  |S12(«)|2).  (2.4) 


Thia  givea  .  .  » 

<  (W)j2.!y  «d)l2.|y  I2.  1 _ •  (2-Sl 

lri:2,“’i  ly2:l'  M  Pi, 21  Su(b»S22(w) 


I  V 


Noting  that  |Cj«*»)IZ  SJU)  and  |G2(W)|2  Su<  <4  art 

ih*  imhi«u  tUnal  power  ipaetn  it  the  rieilvir  In 
Mil  «nJ  IWii,  w  writ#  III*  mafllillwila* 
squared  coherence  ae 


lr 


1.2' 


|Cj(<4!2  !c2«i*!2su(h»2 

[te,  <**  I2  »,«♦*„  «•»]  [fc,c*  f2  suM+ sw2m] 

(1.*) 

Defining  the  eignal-to-noiee  power  in  the  Jth  channel 

»•  . 

|C£U|rS..(«) 


(S/N). 


V"> 


wt  can  writ* 


(S/M), 


(S/M), 


H*)J 


(S.  V 


(S.  •> 


Certain  genern)  etatementa  concerning  thia 
peirwiee  magnitude- .qua red  coherence,  or  just  "co¬ 
herence.  "  cen  new  be  node. 


The  coherence  it  bounded  between  aero  and 

one: 

o  *  lyi,zl2  *  1  •  (».») 

If  the  noiae-to-aignnl  power  goea  to  infinity  in 
either  channel,  the  coherence  will  go  to  aero.  Thia 
will  happen  if  the  aignel  power  in  that  channel  fader 
to  aero.  Alao,  the  noiae-to-aignal  power  ratio  in 
both  channele  moat  go  to  aero  for  the  coherence  to  go 
to  one. 


An  interacting  end  informative  interpretation 
of  the  coherence  between  theee  two  channele  can  be 
made  in  the  following  manner.  Ataume  one  of  the 
channele  ta  noiae  free  (S^lul)  •  0).  Thia  channel 

then  becomea  the  input  eignal.  The  coherence  be¬ 
tween  cnanncla  one  and  two  ia  now  given  by 


, 

!y,  »'  * - ; - - 

ICj Ml  5uM*Syl(U) 


(S/N), 


lefS/NJj  Sj  +  Nj 


(3.10) 


This  magnitude*  a  qua  red  coherence  it  the  fraction  of 
the  power  of  the  output  Xj(t)  which  comet  from  the 
signal  input  patted  through  a  linear  ayetem. 


Since  the  transmitted  eignal  it  generally  not 
available,  it  it  useful  to  look  at  this  physical  inter¬ 
pretation  of  coherence  from  a  different  point  of  view. 
Let  ut  a  imply  take  eignal  *2(t)  ae  our  baeic  eignal  and 
calculate  the  coherence  between  x,(tl  and  our  given 
signal  x~(t).  From  equation  2.7,  our  definition,  thia 
will  be  the  tame  result  ae  if  we  took  signal  x.  (t)  ae 
our  "given"  eignal.  Thue  U(fa>)  in  equation  3.4  ie  re¬ 
placed  by  XgfW): 

Xjfwi  «  H U(«)X2(«)  *  V<12(«)  .  f3.  U) 

ie  the  effective  linear  tranefer  function  be¬ 
tween  output  x^ft)  end  output  Xj(t).  V^j2(t)  ia  lh* 


effective  noise  on  the  transmission  channel.  It  must 
be  noted  that  U.jM  *■  no  longer  nacaaaarily  a  cauaai 
evetem.  N"V  tlta  rnlifrvnri  liatw»an  rliano»la  on* 
and  two  can  be  written  ae  ui  equation  3.  10i 

Ih  (Wlfs-jiw) 

mr - — - r  (s.  i2) 


Irj.jfwtt2  • 


[|H12 


(«)l2S22(tf)+Svtl2(U)'| 


Thia  two-chanti  al  magnitude -a  qua  rad  coherence  ic 
the  retie  of  tha  pewtr  at  output  a,  (t),  which  ia  cauaad 
by  tha  "input  x,(t), "  tranamitt.d  over  tba  affectiva 
linaar  transmission  channel  [H,2(«))  to  Me  total 
power  in  outputs,.  Consideration  of  thia  .ffoctiva 
linaar  tranamiaaion  chanaal  allow,  thia  phyaical  in¬ 
terpretation  of  the  coherence  to  bo  easily  carried 
ovor  to  multiple  coherence. 


4.  Multiple  Coherence  end 
Signal -to-Noiaa  Ratio. 

In  tha  caae  of  M  channel.,  a  ralationahip  batweaa  tha 
input  aignal-to-nolaa  ratio,  and  tha  multiple  coher¬ 
ence,  aimilar  to  the  one  in  the  laat  aection  for  tern 
channel*,  can  be  derived.  The  output  power  epectral 
danaity  cen  again  be  written  in  term*  of  thia  input 
eignal  apactral  density,  the  unknown  channel  tranafar 
functions,  and  tha  effective  channel  noiae  aa 

Su(«)  •  ICjtw)!2  SU(W>  ♦  S¥l(u»  (4. 1) 

ISydall2  -  |GI(U)|2  |GJ(*))|2Su(u)2.  idj  (4.2) 

It  should  be  notad  that  iG^u)]2  S  (<*>)  i*  tha  eignal 
power  danaity  In  the  output  of  theu  ith  channel  and 
S  jfu)  la  the  noise  power  spectral  density  in  tha 
c Kennel.  Tha  general  power  spectral  density  func¬ 
tion  can  then  be  written  aa 

S  (W!  *  S  (w)  C*(W)  GT(e>)  +  Dt*> )  ,  (4.3) 

XX  u 

where 

D(fa>)  ie  e  diagonal  matrix  with  elements  S^.fw), 

end 

GT(*>)  e  [Gj(«)  G2(e>).,.  CM(W)]  .  (4.4) 

'•  * 

The  inverse  required  to  calculate  the  multiple  coher¬ 
ence  from  equation  2. 2  can  now  be  calculated  by 
using  the  following  matrix  inversion  lemma: 

E*  Tl**  -l  .  1  •  T  -l  •  T  -t 
A  +  X  XTj  -A  -A  lX  1  ♦  X^ A  lX  X^A 

(4.5) 

Using  this  lemma,  the  inverse  of  the  spectral 
density  matrix  cen  be  written  (assuming  all  required 
inverses  exist)  as 


i  •  r  t  -i  •  T-i 
(fa))  C  <fa))h+Su<fa)G  (W)D  (w)C  <w)' 

D"'(ai)  . 


S^(U)  «  D*  (*>) 

•  cTi«) 


i 


(4.6) 


*  . 


The  structure  of  this  inverse  cut  be  sees  mere 
clearly  by  noting  that 


Note  that  for  M  equal  te  two  ee  in  section  S  the 
coherence  is  given  by 


t  _i  m  n  lGI(e»)|  9_  fu)  i._  i  i.  i  (S/N)4 

V",G  <w>°  «“>°  <«■».£ -*"5  tj;  tt.7)  lri.*1*  ■>*»»*  -  - 7^2  • 


This  term  is  the  sum  of  ell  output  signal-to-noise 
power  reties.  There/ere,  if,  es  in  the  two. channel 
ceee,  we  define 

(S/N)t  -  |C1(«)|2Sll<«)/Svl(*l)  (4. 9) 

the  inveree  of  the  bracketed  term  in  equation  4. 6  can 
be  written  as 

r  ..i  r  m 

|l  +Su(«)CT(o»D*,(|iI)C*(«)J  el/  1  +£  (S/N)^  (4.4) 

j-1 

With  this,  the  ith  diagonal  element  of  9**(W)  ie 
given  by  5 


-ii,  .  1 

s  («)  - 

vi 


♦  £  IS/Nlj  -  IS/NJj 


M 

14  ^(S/N). 
j-1 


Using  this  and  substituting  equation  4.10  into  equation 
2.  2.  we  find 


7i/1.2,...i.l,!4l,...Ml 


1  -  1  /  Sylu)  S“(w) 


■  a, 

(S/Nlj  > 

l  1 

(S/Nlj  -(S/NJjj 

r  M 

[l  ♦  (S/NlJ 

14  £  (S/N^  -  (S/Nlj 

t  j-1 

Several  special  cases  are  of  interest. 

first  consider  the  situation  in  which  the 
» if nsl •  to-noiee  ratios  in  all  channels  are  the  same: 

(S/Nij  e  (S/N).  «  (S/N)  .  (4.12 

This  fives  the  coherence  of  channel  i  with  respect 
to  the  other  M-l  channels  as 


*1:1,1 _ l-l,  141, , . .  M1 


(S'NriM-1) _ 

[1  -  (S/N)]  [1  ♦  (M-l)  (S/N) ) 


Hew«v«r*  1 1  M  bieemta  v«ry  l»r|«,  th«  coh«r*ac« 

|OM  tO 

*ri:  1 . 2, . . .  i- 1. 141. . . .  M **”“  [I  ♦” (S/N)]  -  M-l*> 

The  formal  requirement  for  this  to  be  valid  ie  for  tits 
signal. to-notee  ratio  and  number  of  channels  te 
satisfy  the  following  Inequality: 

(M  -  1)  (S/N)»  1  .  (4.16) 

However,  based  on  data  from  a  action  3,  equation 
4.  IS  it  identical  te  die  coherence  of  two  channels 
when  one  has  an  infinite  signal -to. noise  ratio  and  the 
other  has  a  signal-to-noite  ratio  (at  frequency  w)  of 
(S/N).  In  this  sense,  a  large  enough  number  of 
weak  channels  (signal. to-noiae  ratio  of  (S/N)]  ie 
equivalent  to  the  sum  of  one  noise. free  channel  and 
one  weak  channel. 

The  second  special  cate  for  equation  4. 11  is 
when  the  ith  channel  has  a  very  large  signal -to-noiee 
ratio.  Letting  (S/N),  become  large  in  equation  4. 11 
and  keeping  all  other  aignal-to-noise  ratios  equal  to 
(S/N)  we  find  that 

|y  p 

,7ltl.w...  .1-1,141,  ...M1 


(M.l)  (S/N) 


(S/N).—  -  [1  +  (M-l)  (S/N)]  ’  ,4-17> 

Note  that  for  low  signal -to-noiae  ratios,  i.e., 

(M-l)  (S/N)  «  1  .  (4.  IS) 

the  coherence  goes  up  linearly  with  the  number  of 
channels  each  is  considered  to  have  the  same 
signal-to-noise  ratio  as  all  others,  i.e.,  (S/N).  As 
M  becomes  larger  or  as 

M(S/N)»  1  .  (4.19) 

this  coherence  goes  to  one  as  it  would  in  the  case  of 
two  noise-tree  channels. 

Next  consider  the  case  where  one  channel  other 
than  tht  ith  channel  has  a  very  high  signal-to-noise 
ratio  relative  to  the  others: 


(S/N)k  >  >  T  (S/N^  -  (S/NIj  -  (S/N)k  .  (4. 20) 


>  . 


I 


1 


\ 


i 

(. 


1  V 


1'nder  these  conditions,  equation  4. 11  1*  ipproxi* 
mxl.iy  Kivan  by 


7i:l,  2, . . .  i-l, i.l 


J»  <S/N)((S/-Nlk 

luiS/N^fuiS/Ng  (4-2U 


The  coherence  between  the  strong  eignel  end  the 
wn*k«r  onua  *»•»•  up  Unfairly  with  dm  mimlmr  of 
w»eb»r  'I'lita  «*»■♦  dm  m< 

th«  M  multiple  coherence  valuer  will  6#  the  ou  in 
which  ihc  Urgeat  ai|MMo*noiM  retie  channel  U 
tired  u  the  reference,  which  i«  *a  expected. 


or 

|yi:1.2...t-l.Wl,...M  |  ~|yt,k  I  .  (4.22) 

.1  it  ell  ocher  channel,  were  not  ueed.  It  (S/N){  ie 
approximately  equal  to  (S/N). ,  thte  meene  that  alt 
weaker  chaaaela  could  be  n.flect.d  and  only  the  two* 
channel  coherence  between  the  (wo  otrongcr  channel, 
could  bo  ueod.  Aho  conotdor  the  case  whan  all 
channela  including  the  ith  have  a  much  lower  signal* 
to-noise  ratio  than  the  kth  channel,  1.  c, , 

(S/N)k  *  '  (S/N)t  .  (S/N)  aU  t  i  k  (4.  IS) 


S.  A  Sample  Statistic  <or 
Multiple  Coherence 

The  true  multiple  coherence  of  o  sot  of  time 
scries  is  c  function  of  the  underlying  statistics  of 
theae  processes.  The  statistics  ore  generally 
unknown  and  must  be  estimated  from  sample  realise, 
liana  of  the  processes.  The  estimates  of  the  baste 
atatictic*  caa  than  be  ueod  te  provide  estimates  of 
die  multiple  coherence  of  the  M  underlying  stochastic 
processes. 


Then,  while  Ihc  coherence  of  the  ith  channel  given  the 
others  is  provided  by  equation  4.  21,  the  coherence 
of  the  kth  channel  given  the  ethers  is 

)v  |2  (M-l)(S/N)k(S/N) 

I  ’k:l,2.  ..k-l,k+l,...Ml  * fiVfS/N) ][1.<M-1)(S/N)J’ 

(4.  24) 

Giving  for  this  cecc 

|rk:  1 . 2 . k*l,k+l,...M  |* 

. -'•'•■-•“I'- 

<4.  25) 

For  the  ceee  of  week  aignele  in  the  other  channels 
tfrorr.  equation  4.  22): 

iyksl. 2, . . .  k-l.k+I , . .  ,m|  ■ 

(4.  26) 

Further  simplify  equation  4.  26  to  the  epacisl  case  of 
(M-IH5/N)  <<1.  (4.27) 

we  have 

|\:1,2 _ k-l,k*l,...M  |  |  (4- 


The  method  of  obtaining  estimaie.  for  true 
multiple  coherence  is  as  follows.  Using  well  docu* 
mented  techniques*’ 4  *  * ",  obtain  sample  eetimatee 
far  each  aleirant  of  the  croeepower  spectral  density 
matrix.  From  chase  sample  estimates 


S  (w)  » 


Sn(N,  Su(-) 


S2,(W)  SH(w» 


S1M<“> 


S2M,“I 


!W“'  SM2,W> 


S  (to)) 
■MM 


<S.l) 


one  calculate*  the  temple  eaumate  for  multiple  co¬ 
herence  in  the  following  manner 

lyi:1.2,...,i-l,ial....M  2*1  -1  /[S^lwis"  <«)].  (5.2) 

where  S^fw)  ie  the  ith  diagonal  element  of  the  inverse 
of  «be  M  -by-M  sample  spectral  density  matrix  S^ftu). 

Det.il.  of  how  to  form  such  eetimatee  are  diacueeed 
a:  length  in  the  literature.  Since  theee  cattnr.atea  ere 
rindom  varieblee  there  has  been  considerable  ttudy  of 
their  distribution.  The  dietributione  of  theee  cross* 
end  ajitopower  spectre!  eetimatee  ere  known  in  closed 
form'. 
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Tha  cla»<l>lorm  expression  for  the  multiple- 

.  .  w  Ilitt 

i.'.c  t->»* u  wi  v<lu««  oi  the  multiple-coherence  test 
luiimt  and  the  ralattva  probability  e<  It*  bala|  In  a 
particular  band.  All  value*  ara  ot  couraa  bounded  by 
zero  and  on*.  The  danalty  lunctlon  la  canditianad  an 
in*  tatal  numbar  oi  dliiarant  time  racorda,  ar  dif¬ 
ferent  atochaatlc  proceaeee,  avallabla  (p).  It  la  alao 
conditioned  on  tha  numbar  ot  Indapandant  aampla* 
available  item  each  of  tha  time  record*  (N).  Thua 
the  denaity  function  of  the  aampla  aatlmata  for  coher¬ 
ence  given  tha  true  coherence  la  given  by* 

P_<y/N.p.|r|2) 

2r,(N.N.p-l<y|yj2) 


(l.Wiy,iN-,p.ii 


Lillie 


p _ (y/N.  p.  |yfS«o  y  >  1  or  y  <0. 


In  aquation  5.  3,  ^Fj(  )  la  the  hypergeometrtc 


Thi*  expreaaion  for  the  denaity  function  of 
multiple  coherence  ta  both  expenatv*  to  calculate  and 
generally  numerically  ill  conditioned.  Thua  to  eval¬ 
uate  the  denaity  numerically,  additional  manipulation* 
are  required.  Creat  difficulty  can  be  encountered  in 
attempting  to  uae  computer  library  expreaaiona  for 
the  hypergeometne  function. 

For  low  value*  of  N  and  p  we  uae  a  tranaform- 
ation  given  in  reference  12: 

2Fl(N,N.p-l:  ly!2  y) 

*  (1  *iylZy)p'1.*2N2Fj(P-i-N.p-li  j  y]2y).  (5.4) 

For  the  caaea  of  intereat,  (p-l-N)  la  a  negative 
integer  ao  that  a  finite  aerite  expanaion  for  thia  latter 
hypergeometne  function  la  available; 

2Flr..N.p.I;|yl2y) 


Other  expreaaiona  for  thia  denaity  valid  to  large  valuee 
of  N  and  apecific  range*  of  true  coherence  |  ^  and  y 
ar*  deecribed  in  reference  1.  There  curve*  of  prob¬ 
ability  of  detection  verau*  probability  oi  fair*  alarm 
are  preaented  aa  ar*  curve*  of  probability  of  detec¬ 
tion  verau*  true  coherence  ter  fixed  value*  of  prob¬ 
ability  of  tale*  alarm. 
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i.'.mg  the**  expreaaiona  for  the  hypergeometric 
:  unction  w*  can  write  for  y  between  aero  and  ones 


